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Abstract 

In a general class of one and two dimensional Hubbard models, we prove upper 
bounds for the two-point correlation functions at finite temperatures for electrons, for 
electron pairs, and for spins. The upper bounds decay exponentially in one dimension, 
and with power laws in two dimensions. The bounds rule out the possibility of the 
corresponding condensation of superconducting electron pairs, and of the corresponding 
magnetic ordering. Our method is general enough to cover other models such as the 
t-J model. 

The Hubbard model and its variants have been attracting considerable interest. But rig- 
orous results are still rare. In one dimension, the Bethe ansatz method has been successfully 
applied [Q both to the ground state and to the finite temperature Gibbs state. In general 
dimensions, Lieb's theorem [Q] and Nagaoka's theorem || on the ground state structures are 
known. In one and two dimensions, Ghosh |3J proved the absence of magnetic ordering at 
finite temperatures. 

In the present letter, we extend McBryan and Spencer's method [Q developed in classical 
spin systems to a general class of Hubbard models in one and two dimensions, and prove 
upper bounds for various correlation functions at finite temperatures. The bounds rule 
out the possibility of magnetic ordering and condensation of electrons or superconducting 
electron pairs such as Cooper pairs or 77-pairs ||. 

We consider a tight binding electron model on the one dimensional lattice Z or the square 
lattice Z 2 M. The Hamiltonian is given by 

H = ~ £ H l ^y c \,c c v,c + V{{n x ,a}) + h * " W 

x,y£Z d <r=Td xeZ d 

with d = 1 or 2. The number operators are defined by n XjCr = cl a c x ^, and spin operators 
by S J X = Ey^.jC^.j'V (-7 = ^2, 3), where (r^ (T ,) (Ji(T / =Ta are Pauli spin matrices and 
cJ. CT , c Xj<J , are the creation and the anihilation operators, respectively, for the electron at site 
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x with spin a. The hermitian hopping matrix (t xy ) is arbitrary, except for the conditions 
that there are finite constants t, R, and \t xy \ < t holds for any x, y, and t xy is vanishing || 
for | re — y\ > R. Note that we can include external magnetic field which is represented by 
complex txy- The interaction V({n XjCr }) is an arbitrary function of the number operators, 
and h x represents local magnetic field or spin-flip impurity. Note that the Hamiltonian (|1|) 
is not necessarily completely isotropic in spin space, but has a global 0(2) symmetry related 
to the spin rotation about the z-axis. We stress that the class of Hamiltonians considered 
here includes not only the well studied models like the (standard) Hubbard model or the 
periodic Anderson model, but also many of their variants with , e.g., long-range, random or 
spin-dependent interactions. 

To define the Gibbs state, we replace the infinite lattice with a finite lattice of linear di- 
mension L with periodic boundary conditions. The thermal expectation value of an arbitrary 
operator A is defined by 

(A) Tl{Ae ~ 0H) (2) 

where the trace is over all the electron states. We consider the infinite volume state defined 
by 

(A) = lira (A) L , (3) 

with the electron density fixed to p. Our result is independent of p and thus applies to grand 
canonical averages as well. 

The main result of the present letter is the following. 

Theorem: There exist finite constants || a, 7, 5 and a function /(/?) which depend only on 
the hopping matrix (t xy ). The function is decreasing and behaves as /(/3) w for 
(3 ^> 5 and f(/3) ~ (2/5) \ ln/3| for (3 < 5. In a two dimensional model in the class described 
above, we have 

i T ii c ^T%l + H.c.) < 2\x-y\-«M\ (4) 
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4^ + H.c.) < 2\x-y\' a ^\ (5 ) 

for any finite f3 and for any x, y with sufficiently large \x — y\. If the local field has the form 
lXj. = (0, 0, h x ) we further have 

^lSy + S^Sy^ < \x — y\ a f(@\ (6) 

for any finite (3 and for any x,y with sufficiently large \x — y\. In a one dimensional 
model, we have the above bounds (f|), (0) and (^) with the right-hand-sides replaced with 
2exp[-7/(/3)|z -y\], 2 exp[->y f (2(3)\x -y\/2] and exp[-jf((3)\x - y\), respectively. 
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The above bounds rigorously rule out the possibility of the corresponding condensations 
of electrons or electron pairs and of the corresponding magnetic ordering. The bound (§), 
for example, inhibits the condensation of singlet electron pairs such as the Cooper pairs 
or the ??-pairs ||. However our method can be easily extended to rule out any kind of 
condensation which is related to a spontaneous breakdown of the quantum mechanical global 
£7(1) symmetry. It is also straightforward to extend the method to cover other systems such 
as the Hubbard model with nonlocal spin-flip term or the t-J model ]L(J. The explicit upper 
bounds for the correlation functions provide further information about the propagation of 
electrons, electron pairs and magnons. The astonishing generality of the theorem, especially 
the complete arbitrariness of interactions, may be regarded as a sharp demonstration of the 
fact that the electron hopping plays a fundamental role in various condensation phenomena 
in itinerant electron systems. 

The power law decaying upper bounds in the theorem are certainly not optimal at high 
temperatures, where one generally expects to have exponential decay. Even in low tem- 
peratures, a class of models which are sufficiently close to the antiferromagnetic Heisenberg 
model is expected to show exponential decay. Among the varieties of models covered by 
the theorem, however, one might well find those which exhibit "exotic" phase transitions 
leading to power law decay. It is notable that the power indices in the upper bounds (f|), 
(§), (H) are proportional to (3~ l at low temperatures. This means that the slowest possible 
decay in these models is of the Kosterlitz-Thouless type. In one dimension, the exponentially 
decaying upper bounds in the theorem provide upper bounds for various correlation lengths. 
The bounds, which are proportional to (3~ l at low temperatures and to | ln/3| at high tem- 
peratures, reproduce a typical crossover behavior of correlation lengths in one dimensional 
tight-binding electron systems. 

Our proof is based on the method developed by McBryan and Spencer |J for classical 
spin systems, and on its extension to quantum spin systems by Ito ||11||. In these works, the 



global continuous symmetry of the spin space played an essential role Jl2[ . Our strategy here 
is to make use of the global U(l) symmetry related to the quantum mechanical phase. In this 
approach, we do not have to make further assumptions on the symmetry of the system since 
the U(l) symmetry exists in any quantum particle systems. We believe that the present 
method can be extended to much larger class of quantum particle systems. In the present 
letter, we restrict ourselves to the lattice fermion problems, which are free from ultraviolet 
divergence. 

The absence of magnetic ordering in one and two dimensions was proved by Ghosh @], 
who extended the Bogoliubov inequality method of Mermin and Wagner's We note 



that, by combining the Mermin- Wagner argument with the idea to make use of the quantum 
mechanical U(l) symmetry, one can also prove the absence of condensation of electron pairs 
(or electrons). To do this, one should replace the operators A and B in || with the Fourier 
transforms of the number operator n x = n x ^ + n Xi i and of the order variable O x = c X) ^c x ^ 
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(or c Xj0 ), respectively. We also note that the Mermin- Wagner argument can be extended to 
cover non-translation-invariant models as those considered here. 

In what follows, we describe the proof of the bound (f|) in detail. We first prove the 
bound in a finite periodic lattice of linear dimension L, and then take the limit L — > oo. 
To make use of the global quantum mechanical symmetry, we note that the U(l) gauge 
transformation is represented by the unitary operator 

G(6) = Hexpi-i0 u n Uta ], (7) 

U,(T 

where 9 = {9 U } is an arbitrary real function on the lattice. In the following, however, we let 
9 U to be pure imaginary, in which case the operator G{9) is no longer unitary. Since G(9) is 
invertible, we have 

Tr[Ae-P H ] = Tr [G{9)AG{9)- 1 exp[-/3G(9)HG(9)- 1 ] 

for arbitrary complex 9 U . Here the transformed Hamiltonian is 

G(9)HG(9)- 1 = -J2 tu, v e-^-^4 jCT cv + V({n U;a }) + £ h « • S«- (9) 

U,V,IJ u 

Let ip = {ip u } be a real function which will be specified later. We consider the operator 
G(—i<p) obtained by setting 9 = —iip = {— i<p u } in ([?]). Let us fix lattice sites x, y, and take 
A = c i,T c i,l c y:t c s/,l- Straightforward calculations show 

G (-tip) AG (-tip)- 1 = exp[-2(^ - <p v )]A, (10) 



and 
where 

and 



G(-ip)HG(-ip)- 1 = H + U + iP, (11) 
U = - £ t u , v {cosh(<p u - <p v ) - l}cl a c v ^, (12) 



P = -2 £ tu,v sinh(p u - </?„)4 )CT c„ j0 . (13) 

u,v,a 

are hermitian matrices. 

We can bound the right-hand-side of (§) as 

Tr [G(-z^)AG(-'^)- 1 eyv[-f3G(-ip)HG(-ip)- 1 } 

< e -2(^-^)Q|^* y4 || oo ^ ) l/2 Tr | e -(/3G'HG- 1 )/2 e -(/3G- 1 HG)/2j 

< e -2(^-^) Tr [ e -/3(H+(7)j 

< e-^-^^e^HooTrfe^], (14) 
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where ||0||oo denotes the maximum of the absolute values of the eigenvalues of a hermi- 
tian matrix 0. To prove the above bounds, we use the following inequalities for oper- 
ators (matrices) on a finite dimensional space, i) The Schwartz inequality; Tr [OP] < 
{Ti[0*0]Tr[P*P]} 1/2 with 0,P arbitrary, ii) |Tr[OP]| < ||0||ooTr[P] with O hermitian 
and P positive, iii) Ti[(0*) N N ] < Ti[(0*0) N ] withiV = 2 m and O arbitrary (|TJ. iv) The 
Golden- Symanzik-Thompson inequality |15] Tr[e° +P ] < Tr[e e p ] where 0,P hermitian. To 
show the first bound in (0), we set W = exp[-/3G(-tip)HG(-tip)- 1 /2], A' = G(-i V )AG(-i V )-\ 
and use i) and ii) to get 

Tr[(A'W0iy] < (Ti[A'*A f WW*]Tr[WW*]) 1/2 

< (WA^A'W^Y^TrlWW*]. (15) 

The second bound follows by noting that ||^4||oo = 1, and setting 

(3G{-i V )HG{-i V )- v 



X = exp 



2N 



(16) 



UV I ; 



to get Tr[X (X*) ] < Tr[(XX*) ] from iii). The right hand side converges to 
Tr[exp[— (3(H + U)}} as N — > oo. The third bound is an easy consequence of iv) and ii). 

Now we choose (p. Let X uv be real hopping matrix elements that satisfy X uv = X vu > \t, 
and \ uv = for \u — v\ > R. We further require X uv to be periodic, i.e., there are positive 
integers p, q, and X uv = X u+ d !V +d holds for any d = mpei + mge 2 where m, n are arbitrary 
integers and ei,e 2 are two unit vectors of the lattice. (In one dimension, we of course set 
d = mpei.) We assume that the lattice size L is a common multiple of the periods p, q. The 
conditions imposed on t uv ensures the existence of such X uv . (The simplest choice, which is 
always possible, is X uv = t for \u — v\ < R, and X uv = otherwise. By choosing X uv which is 
"closer" to t uv , however, one gets better constants in the resulting bounds.) Let / = {/«} be 
a function of the lattice sites, and define a lattice Laplacian A by (A/% = J2 V ^uv(fv — /«)■ 
We let <f = {f u } be the unique solution |L6] of the Poisson equation — (A<p) u = q(S XtU — 5 y}U ) 
with a zero-point condition ip y = 0. The "charge" q > will be determined later. By using 
the periodicity of X uv and explicitly writing down the solution in terms of the Fourier series, 
one finds that (p has the following two properties |J. PI) There exists a finite constant S, 
and \tp u — tp v \ < q5 holds for any u, v with \u — v\ < R. P2) In the L — > oo limit, one has 
fx > Ql\ x ~ y\ in one dimension and ip x > gain \x — y\ in two dimensions for sufficiently 
large \x — y\ with finite constants 7, a. 

Noting that the above property PI) implies cosh(yj u — (p v ) — 1 < g(q)(<p u — tp v ) 2 for 
\u — v I < R with g(q) = {cosh(g5) — l}/(qS) 2 , we have 



exph/^IU < 



P \t U v[cosh((p u - (f v ) - 1] 
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< exp 



exp 



-2[3g(q)J2<Pu(^)u 



where we have used ||<4 ^.c^ + ci„c 



V (j u*& II oo 



exp {2f3g(q)qip x } , (17) 
1 to get the first bound. By substituting the 



bounds (fH|) and ( |T7| ) into we get 

Tr[Ae-P H ] 



Tr[Ae-P H ] 
To optimize this bound, we define 



< exp[-2(p x + 2(3g(q)q(p 0l 



f((3) = max[2g - 2/?cT 2 {cosh(g<5) - 1}], 

q>0 



(18) 



(19) 



which is manifestly decreasing in /3, and has the asymptotic behavior stated in the theorem. 
By using the property P2) of ip and letting the "charge" q to be the maximizer in the above, 
we finally get 



< 



f exp[-/(^)7|ar-y|] 
lexp[-/(/3)aln|x - y\] 



\x - y\ 



if d = 1; 
-am if d = 2, 



(20) 



for sufficiently large \x — y\. Thus the bound (^) has been proved. 

The bound ([5]) is proved in exactly the same manner. To prove the bound (y), we set 
A = S^Sy, and perform a spin-dependent unitary transformation represented by G{6) = 
Y[ u ,cr exp[— ia9 u n UiCT }. The rest of the proof proceeds in exactly the same way as the above. 

We wish to thank Elliott Lieb for letting us know the reference Q. 

Note added (September 1997): In |T7|, extensions of our results to models with long range 
hopping was discussed. 
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